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Spin and thermal conductivity in classical disordered spin chain 
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Transport quantities of the classical spin chain with the quenched disorder in the antiferromagnetic coupling 
,Ji are evaluated using the dynamical simulation at finite temperatures T > 0 . Since the classical model is 
nonintegrable, spin and thermal conductivities remain finite even in the pure case. On the other hand, the role 
of disorder becomes crucial at low T leading to a vanishing transport due to the Anderson localization within 
the linearized regime. The crossover from the insulator to the conductor appears both for the spin and thermal 
transport at quite low T* J. Still the many-body localization regime at T > 0 evidenced by extremely short 
mean free paths can be strongly enhanced by introducing into the model an additional staggered field. 
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I. INTRODUCTION 

Properties of low-dimensional spin systems have been in¬ 
vestigated intensively since many decades. Since the one¬ 
dimensional (ID) S' = 1/2 antiferromagnetic (AFM) Heisen¬ 
berg model is the first quantum many-body model consid¬ 
ered theoretically and its ground state properties are exactly 
solvable by the Bethe Ansat^l it has been the playground 
for numerous analytical and numerical approaches. In recent 
decades, experimental realisations of the spin-chain physics 
has been found in several classes of materials, the best and 
nearly ideal example of the isotropic Heisenberg S' = 1/2 
model being cuprateA particular signature of the ID 
S = 1/2 model and its integrability are transport properties 
which are anomalouJ^. It has been shown that the heat trans¬ 
port is dissipationless at any temperature T > (J^^, at the same 
time also the spin conductivity (diffusivity) is infinite (ballis¬ 
tic) in the anisotropic easy-plane regime (anisotropy A < 1) 
at any T > cP. 

In contrast to clean spin chains the transport in disordered 
systems is much less understood and represents a challenge 
and very active theoretical topic. It is well established that in 
a quantum S' = 1/2 AFM Heisenberg chain and more gen¬ 
erally in the system of ID repulsive fermions at T = 0 any 
random quenched disorder induces localisation and the van¬ 
ishing of transport coefficient^® This is the consequence of 
the Anderson localisation phenomenorP^ti^ in ID persisting 
or even enhanced in the presence of interaction within the 
ground state. The scenario of the many-body (quantum) lo¬ 
calisation (MBL^P^JG^which would manifest itself with the lo¬ 
calisation and no d.c. transport at T > 0 is investigated inten¬ 
sively at present, mostly within the random Heisenberg model. 
So far the evidence for the latter has been found in numerical 
solutions of strongl y dis ordered S' = 1/2 AFM model with 
random local fields^IE], while some studies indi cate more 
on a crossover into a nearly localised regimJSC^. We note 
that there are recently experimental realisations of disordered 
(but isotropic) spin S = 1/2 chains BaCu 2 (Sii^Gea;) 07 n^, 
Cu(py) 2 (Cli_a:Bra;) 2 ® and Ca-doped Sr 2 CuoJ^ where the 
disorder enters predominantly via random exchange couplings 
Ji and has been theoretically studied in connection with the 
concept of random singletPBlSI 


The classical AFM spin-chain model, which can be re¬ 
garded as the S —t oo limit of the quantum model, is nu¬ 
merically much easier to deal with. Without disorder the spin 
dynamics of is nonlinear and nonintegrable, hence one obtains 
at T > 0 finite transport coefficient^, unlike to S = 1/2 
modeP but on the other hand closer to S' > 1 case^. Still, 
the role of disorder bears some analogies to the MBL when 
considering low but finite T > 0. Namely, in the low T 
regime one can imagine the transport being dominated by lin¬ 
ear excitations whic h exhibit the Anderson localization in a 
random systenJAliSl Numerical studies of the classical model 
so far performed at T —> ocPdid not reveal the absence of 
diffusion. In the following we perform the numerical study 
of the d.c. spin conductivity as and the thermal conduc¬ 
tivity K within classical spin model in the whole T regime. 
We show that indeed at low T —> 0 both coefficients van¬ 
ish due to the localization phenomena. Here the fundamental 
question is whether the (classical analogue) of the MBL per¬ 
sists in a finite T window. We find the evidence only for a 
insulator-conductor crossover regime, in this sense consistent 
with studies of transport in disordered anharmonic chain^. 
Still, the quasi-insulating regime characterized by extremely 
short spin and thermal mean free paths (MFP) ^ 1 and 
the crossover temperature T* can be strongly extended by in¬ 
troducing a constant staggered field. In spite of strong T de¬ 
pendence, both transport MFP qualitatively satisfy the simple 
Wiedemann-Franz law® 4 ^ It- 


The paper is organised as follows. In Sec. II we introduce 
the disordered classical ID spin-chain model. We define the 
appropriate spin and thermal currents as well as the linear- 
respond equations for the d.c. spin conductivity as and the 
thermal conductivity k. In Sec. Ill we describe the numeri¬ 
cal protocol for the the initial spin configuration correspond¬ 
ing to the equilibrium at chosen temperature T and further 
simulation of the spin dynamics and calculation of d.c. trans¬ 
port properties. Results for as{T) and k(T) are presented in 
Sec. IV for various disorder and the whole range of T. A 
particular analysis in terms of corresponding mean free paths 
Is, h, respectively, is given in Sec. V and discussed in relation 
with the Wiedemann-Franz ratio W = It/h- 
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II. MODEL AND TRANSPORT QUANTITIES 


III. NUMERICAL PROCEDURE 


We consider in the following the ID AFM isotropic model 
of classical spins (rotors) = S'e^, 

= -S,, (1) 

i i 


where are unit vectors and (exchange) couplings Ji > Q 
are antiferromagnetic on all bonds but random, in analogy 
with the S — \/2 casd^^®l We also choose Ji randomly 
distributed within an interval J — dJ < Ji < J + 6J (with 
SJ < J so there is no possibility of cutting a chain) and are 
as well uncorrelated between sites. Since we fix furtheron 
S' = 1, J = 1 and the lattice parameter a = 1 the only rel¬ 
evant parameter (in the absense of local fields, Bi = 0) is 
the disorder strength 0 < (5 J < 1. Our study deals with the 
disorder in the off-diagonal Ji, the motivation being also in 
the material realisations where random fields are hard to 
justify. Nevetheless, we also comment on the influence of the 
staggered field B^ = (—l)*i3ez which can, e.g., emerge as a 
mean field due to interchain couplinj^. We show that B 0 
has a very strong effect on transport at T —0. It should be 
also reminded that the thermal transport in the classical model 
with random B^ has been studied for T ocP. 

The equations of motion for the classical model, Eq. Q. 
are 

JO <9// 

^ = -S.X — = -S,x(J,_iS,_i + J,S,+i-B,). (2) 
at ooi 


Our aim is to study the spin current jg and the e nergy current 
jE which are defined via the continuity equation d'^ l ^^ l. 


dSf 


= -07+1 -jT)> 


dt 

where hi is the local energy 
= 2 ['^i —1 ■ 


dh, 

dt 


= -0^i-jf), 


(3) 


. JA ■ S,+i] - B, • S,, (4) 


The explicit expression for both currents then follow, 

AS ^ J( CX qy _ qy qx \ 

Ji — 5 

jf = i J._i[J,_2(S,_2 X S._i) • S, -f J,(S,_1 X so • S,+i] 


- hJi-1 (Si_i X Si) • (Bi_i 


B, 


(5) 


The d.c. spin conductivity ag and the thermal conductivity k 
are within the (Kubo) linear response given through correla¬ 
tion functions^, 

T 

= J {J%t)J%0))dt, 

0 

r 

K=j^lim J {j^{t)J^{0))dt, 


( 6 ) 


where J'* = 

ergy currents, respectively. We also put fcs = 1 so that 
T = TkB/JS"^, (Jg = ag/[Sa?) and k = n/lksdSa^) are 
all dimensionless. 


Thermal values of agiT) and k{T) for given T and disor¬ 
der 5J are calculated numerically in two steps. First, we use 
modified Landau-Lifshitz-Gilbert (LEG) equations for spin 
dynamicJS] in order to thermalise the system to given T, i.e., 
to find the T-dependent initial spin configurationPH, 


d^ 

dt 


1 + a'^ 


X 


6 


^-aS.X 14 


dH \ 


( 1 ) 


with a damping parameter a, and ^i representing random vec¬ 
tors due to r > 0 thermostat, with properties of the Gaussian 
white noise. 


(6(i)) = 0, (6(ii) • e,(t2)) = 2aT5,,5[h - ^ 2 ), (8) 

In the concrete numerical realisation we perform the cal¬ 
culation on L = 10000 sites. The LEG thermalization, 
Eqs. (|7]l,([^, we independently check via Boltzmann thermo¬ 
dynamic value of the energy ^^(r) = {H). Clearly for T <C J 
longer thermalization times tg are required so that we employ 
fg 1000. Furtheron final LEG spin configurations are used 
as initial conditions for dynamical spin evolution via Eq. Q) 
where we use the standard fourth-order Runge-Kutta with the 
small time step St ^ 0.01 (tested by the conservation of spin 
norm S' = 1 and the total energy E) and the evolution times 
up to T ~ 10000. Kubo formulas, Eqs. (|^, are then used to 
evaluate <Jg(T) and k(T) at various disorders SJ. 

Such a protocol for the evaluation of transport properties, 
standard for generic (ergodic) systems, should be carefully re¬ 
considered and tested for our model, being a candidate for the 
MBL phenomenon. In the case of MBL one could expect in 
an isolated system the lack of thermalisation, dependence of 
the (long-time) linear response on the initial conditions etc. 
Still, as shown later our results do not confirm the strict MBL 
behavior (except at very low T) and consistently the initial 
conditions do not influence the final result. Still the closeness 
to MBL manifest itself in long thermalisation and evolution 
times to, r ^ 1. 


IV. RESULTS: SPIN AND THERMAL CONDUCTIVITIES 

In Fig. 1 we present numerical results for the T dependence 
of the renormalized dg = Tag and k = T^k. It follows from 
Eqs. (j^ that such redefined dg, k remain finite (and constant) 
at T —>^ 00 . Obtained results already reveal general features 
of the d.c. transport in pure and disordered classical chains: 

a) In contrast to the quantum S' = 1/2 spin-chain modeP in 
the classical model both ag{T) and n{T) are finite for T > 0, 
even without disorder at SJ = CP^. 

b) From Fig. 1 it is evident that d{T) and k{T) are both ap¬ 
proaching constants at T > 0.5 for any disorder SJ. Our re¬ 
sults for the clean case SJ = 0 for kao = k{T —>■ 00 ) ^ 0.55 
are consistent with the one previously obtainecP^ 

c) Instead of spin conductivity ag more often considered quan¬ 

tity is the s^ diffusivity Ds,+ determined by the relation 
Dg = (Ts/xJ^where /LT is the uniform spin 









3 



<5J=C 

J 

.8, T=0.02 

L 

1 

F 


2000 2400 2800 

i 


Figure 2. (Color online) An example of local spin deviations emerg¬ 
ing from the linearized Eqs. <2} in case of 5J — 0.8 and T ~ 0.02. 


V. TRANSPORT MEAN FREE PATHS AND 
WIEDEMANN-FRANZ RATIO 


Figure 1. (Color online) Temperature dependence of renormalized: 
a) spin conductivity Bs = Tcts, and b) thermal conductivity k = 
T^k for various disorders 5J — Q — 0.8. 


susceptibility. Taking into account Xs(T —>■ oo) = 1/(3T) 
and our value from Fig. la — ^s{T —>■ oo) ~ 0.25 we 
would get the limit Dg{T —>■ oo) ^ 0.75 which can be well 
compared with the high-T expansion result for the classical 
case Ds{T ^oo) = (2/5)v^l0/3 = 0.7pl (note the correc- 
tion in Ref.®). 

d) It is plausible that the disorder SJ>0 reduces both koo and 
(Too. While for weaker disorder, e.g. SJ = 0.2, the T ^ oo 
transport is not much effected by the randomness, the depen¬ 
dence is more pronounced for SJ > 0.4, in particular for kao- 
The latter difference can be partly traced back to the explicit 
forms of and jf, Eqs. S, where the exchange coupling 
and their disorder enters linearly in j®, but quadratically in 


f) Both <Js{T) and k,{T) are not diverging or approaching fi¬ 
nite values, but rather vanishing in the regime T ^ Q (more 
evident in later plots). The latter T window is quite narrow 
but numerically well resolved. The absence of transport at 
T —>^ 0 is a direct indication of the localization phenomenon. 
Namely, Eqs. Q can be linearized around the AEM ground 
state S° = (—l)*e^. Resulting eigen-solutions of linear equa¬ 
tions, Eq. 0, for Sf (or Sj) are then localised according to 
the Anderson localisatiorinnillll jn ID systems. To confirm 
the onset of localisation we present in Eig. 2 a typical exam¬ 
ple of local spin deviations 5S^ obtained as the solution of 
the linearized Eqs. for large disorder SJ = 0.8 and low 
T = 0.02. It is, however, evident (see later) that the validity 
and feasibility of such a linear approximation is restricted to 
very low T ^ 1 at B = 0, but can be substantially enhanced 
with the staggered filed B > 0. 


It is instructive to analyse the transport data (with the em¬ 
phasis on low T J ) through the phenomenology of 
the standard kinetic theory. One can introduce the concept 
of thermal (heat) MEP It which is in our case the transport 
MEP of AEM magnons representing the relevant low-T ex¬ 
citations with a constant magnon velocity v = 2J (emerg¬ 
ing form the classical AEM dispersion Cq = 2Jsing). It 
then enters the thermal conductivity as k = Cyvlt where 
Cv = {l/L)dE/dT is the specific heat of the system. On 
the other hand, the spin-diffusion MEP Ig can be most reason¬ 
ably extracted via the diffusion constant Dg = vlg- In general 
both MEP are different ls{T) ^ kiT). 



T 

Figure 3. (Color online) Specific heat Cv and uniform susceptibil¬ 
ity Xa VS. T for the pure case (5J = 0 and different values of the 
staggered field B. 

In Fig. 3 we first present Cv{T) as well as Xs{T), needed to 
evaluate It and Ig, respectively. Their dependence on disorder 
SJ (not presented) is quite weak. More significant and impor- 
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tant for further discussion is the dependence on the staggered 
field B. We note that in contrast to the well known S' = 1/2 
quantum case at B = 0 the classical model has a different 
low-T dependence for Cy ^ ks at T J, while for T > J 
one gets Cy ~ 1/(3T^). On the other hand, the behavior of 
Xs{T) has at i? = 0 a similarity to the ID S = 1/2 AFM, 
approaching a constant for T —>^0. 
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T 


Figure 4. (Color online) Transport mean free paths vs. T: a) spin 
diffusion (s(T), b) the thermal and c) the Wiedemann-Franz 

ratio W = It/Is for different disorder strength SJ = 0 — 0.8 and 
B — 0. Note the log scale. 

In Fig. 4 we show (note the log scale) the extracted T de¬ 
pendence of the spin and thermal MFP 4 respectively, for 
B = 0 but different disorder 0 > 6J < 1. The qualita¬ 
tive and even quantitative behavior of both quantities is quite 
similar which is not surprising since the equality h ^ Is 
would indicate the validity of the Wiedemann-Franz lawEH 
even for this system. First we note that for a pure system 
SJ = 0 both MFP increase and diverge by lowering T —0 
which is a normal behavioil^ emerging from restricted scat¬ 
tering processes at T —>^0. At higher T > 0.1 finite dis¬ 
order SJ > 0 reduces both MFP, but does not change the 
qualitative behavior. In contrast, at low T < T* whereby 
the crossover T* we define by = max. the Ander¬ 

son localization mechanism sets in and the trend turns, i.e., 
h,t(T —>■ 0) —?► 0. It is plausible that T* is increasing with 
SJ, still it is quite puzzling that the crossover appears very 
low, typically at T* < 0.05 ^ 1. This implies that even at 
T > 0.1 nonlinear excitations become dominant over linear 
excitations, having analogy to AFM magnons. It should be 


also noted that our results reveal (at B = G) even at mod¬ 
est SJ > 0.4 a narrow but numerically well resolved regime 
T < 0.01 < T*, where both MFP become smaller than the 
lattice spacing ^ 1- 

In order to enhance the localization regime, we consider 
further the influence of the staggered magnetic field B ^ 0. 
Such a field has a great impact on chain properties, e.g., in the 
quantum case this would lead to a spin gap in the excitation 
spectrum. As evident from Fig. 3 i? 7 ^ 0 has an to some 
extent analogous effect also in the classical chain, in particular 
Xs{T —>■ 0) —^ 0 as well well as Cy{T —^ 0) —^ 2 for S 7 ^ 0. 
But the most dramatic effect is on the transport properties. In 
Fig. 5 we present results for ls,t{T) obtained for B = 0 — 
0.5 at fixed disorder SJ = 0.4. It is evident that B > 0 
enhances the localization crossover temperature T* to values 
T* cx B. Moreover, the values of 4 t are strongly reduced 
at low T < T* even at modest SJ and small B = 0.1 where 
the maximum at T = T* hardly goes beyond the ’minimum 
metallic one’, i.e. Ig^t > 1. In B > 0.1 the transport MFP are 
clearly below the naive localization criterion at any T > 0. 
The origin of the localization regime enhanced by B > 0 
is not hard to rationalize, since the introduction of B 7 ^ 0 
extends the validity of the linear approximation to Eqs. to 
much higher T cx B and consequently also to the enhanced 
manifestation of the Anderson localization. 
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Figure 5. (Color online) Transport mean free paths vs. T: a) spin 
diffusion ls{T), b) thermal lt{T), and c) the ratio W = h/h for 
different staggered fields B = 0 — 0.5 and fixed disorder SJ = 0.4 

In connection to results in Figs. 4,5 it is also interesting 
to comment on the relation between both MFP: ls(T) and 
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lt{T). Well known Wiedemann-Franz lavJ^, mostly valid for 
transport in metals, translated to our system would require an 
equality of both MFP It ^ 1^ in the low (i.e. T < J) regime. 
Indeed we find that in the whole range of T presented in Fig. 5 
the ratio W — U/ls ranges 0.5 < W < 2 in the whole pre¬ 
sented T regime, and this in spite of wide range of actual val¬ 
ues of 10“^ < /g t < 10. This implies that the concept of 
constant W makes sense not just at modest scattering but even 
in the quasi-localized (nearly-localized) regime. 

VI. CONCLUSIONS 

In conclusion, presented numerical results for the spin and 
thermal transport in classical disordered spin chain reveal 
some similarities as well as differences to the transport in 
quantum spin chains. First of all, numerical studies for classi¬ 
cal spin models can be performed for large systems reaching 
reliable results for transport quantities even for low T J 
and arbitrary disorder which remains a severe challenge for 
quantum S' = 1/2 model^, even more for larger 1/2 < S < 
oo. 

In the pure case 5J = 0 the S = 1/2 model is spe¬ 
cific due to the integrability and consequently the ballistic 
transport. He nce o nly the disorder can serve as the scatter¬ 
ing mechanisirP^EIl On the other hand, in the classical model 
without disorder <5 J = 0 the mean-free paths remain finite for 
any T > 0, with diverging Lt, Is ^ oo for T 0, and in this 
sense qualitatively approaching the dissipationless transport 
in pure S = 1/2 modelJ^. Still it should be reminded that 
there are qualitative difference between S = 1/2 and classical 


spin model at T —> 0 in static quantities, in particular in the 
specific heat C'v(T). 

The most interesting common aspect of ID classical and 
quantum spin S — 1/2 transport is the onset of localiza¬ 
tion due to disorder SJ > 0 in both models at T —>^0. It 
is evident from our results, that there exists for any SJ > 0 
a characteristic T*, where for T < T* transport MFP It, Is 
decrease and eventually vanish with T —> 0, which can be un¬ 
derstood with the dominant role of the linearized equations of 
motion, Eqs and the Anderson localizastion of the eigen- 
solution. On the other hand, for T > T* the transport is 
already dominated by nonlinearity effects. Our study shows 
that within a classical model there is no abrupt transition be¬ 
tween both regimes but rather a crossover at T ^ T* where 
mean-free paths ( are both maximum. It is remarkable that 
in a isotropic model the crossover, which indicates the onset 
of nonlinear excitations, appears already at very low T* <C J. 
Such a low-T regime could be hardly reached in numerical 
analysis of quantum models. On the other hand, our results 
confirm the finding in the quantum S — 1/2 case that the 
introduction of local fields Bi, both staggered as well as the 
random ones, can induce even stronger effect on the trans¬ 
port propertie d' This brings our model and results closer 
the presumable many-body localizatiorfSEUbut still as an ap¬ 
proximate description. 

It should be also noted that our results at lower T < J agree 
with the qualitative validity of the Wiedemann-Franz law, i.e. 
equality of MFP It ~ Is, recovered even in nearly localised 
regime Is^t ^ 1- 

This work has been supported by the Programs P1 -0044 and 
Pl-0112 of the Slovenian Research Agency (ARRS). 


‘ H. Bethe, Z. Physik 71,205 (1931). 

^ A. V. Sologubenko, K. Gianno, H. R. Ott, U. Ammerahl, and A. 

Revcolevschi, Phys. Rev. Lett. 84, 2714 (2000). 

^ N. Hlubek et al, J. Stat. Mech. 3, 3006 (2012). 

* for a theoretical review see, e.g., X. Zotos and P. Prelovsek, in 
Strong Interactions in Low Dimensions (Kluwer Academic Pub¬ 
lishers, Dordrecht), p.347 (2004). 

^ X. Zotos, F. Naef, and P. Prelovsek, Phys. Rev. B 55, 11029 
(1997). 

® A. Kliimper and K. Sakai, J. Phys. A 35, 2173 (2002). 

^ X. Zotos, Phys. Rev. Lett. 82, 1964 (1999). 

* C. A. Doty and D. S. Fisher, Phys. Rev. B 45, 2167 (1992). 

^ P. Schmitteckert, T. Schulze, C. Schuster, P. Schwab, and U. Eck- 
em, Phys. Rev. Lett. 80, 560 (1998). 

P. Anderson, Phys. Rev. 109, 1492 (1958). 

" P. Dean, Proc. Phys. Soc. 84, 727 (1964). 

for a review see B. Kramer and A. MacKinnon, Rep. Prog. Phys 
56, 1469 (1993). 

D. M. Basko, 1. L. Aleiner, and B. L. Altshuler, Ann. Phys. 321, 
1126 (2006). 

V. Oganesyan andD. A. Huse, Phys. Rev. B 75, 155111 (2007). 

J. H. Bardarson, F. Pollmann, and J. E. Moore, Phys. Rev. Lett. 
109, 017202 (2012). 

A. Karahalios, A. Metavitsiadis, X. Zotos, A. Gorczyca, and P. 
Prelovsek, Phys. Rev. B 79, 024425 (2009). 


O. S. Barisic and P. Prelovsek, Phys. Rev. B 82, 161106(R) 

( 2010 ). 

** 1. Tsukada, et al. Phys. Rev. B 60, 6601 (1999). 

M. Thede et al. Phys. Rev. B 86, 180407 (2012). 

“ A. Mohan et al, Phys. Rev. B 89, 104302 (2014). 

S.-k. Ma, C. Dasgupta, and C.-k. Hu, Phys. Rev. Lett. 43, 1434 
(1979). 

D. S. Fisher, Phys. Rev. B 50, 3799 (1994). 

J. Herbrych, J. Kokalj and P. Prelovsek, Phys. Rev. Lett. Ill, 
147203 (2013). 

J. Kokalj, J. Herbrych, A. Zheludev, and P. Prelovsek, Phys. Rev. 
B 91 155147 (2015). 

A. V. Savin, G. P. Tsironis, and X. Zotos, Phys. Rev. B 72, 140402 
(2005). 

J. Karadamoglou and X. Zotos, Phys. Rev. Lett. 93, 177203 
(2004). 

V. Oganesyan, A. Pal, and D. A. Huse, Phys. Rev. B 80, 115104 
(2009). 

A. Dhar and J. L. Lebowitz, Phys. Rev. Lett. 100, 134301 (2008). 

N. W. Ashcroft and N. D. Wagner, Solid State Physics, (Holt - 
Saunders Int. Ed., New York, 1976). 

R. Kubo, M. Toda and N. Hashitume, Statistical Physics, Springer 
Series in Solid- State Science Vol. 2 Springer-Verlag, Berlin, 
(1991). 

" J. L. Garcia-Palacios and E. J. Lazaro, Phys. Rev. B 58, 14937 







6 


(1998). ” D. L. Huber, J. S. Semura, and C. G. Windsor, Phys. Rev. 186, 

D. Forster, Hydrodynamic Fluctuations, Broken Symmetry, and 534 (1969). 

Correlation functions, (Frontier of Physics, ed. D. Pines) (1975). 

“ D. L. Huber, Prog. Theor. Phys. 39, 1170 (1968). 



